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SPINS OF THE SUPERMASSIVE BLACK HOLE IN M87: 

OBSERVATIONS 
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ABSTRACT 

The rapid TeV 7-ray variability detected in the well-known nearby radio galaxy M87 implies an extremely 
compact emission region (5-10 Schwarzschild radii) near the horizon of the supermassive black hole in the 
galactic center TeV photons are affected by dilution due to interaction with the radiation field of the advection- 
dominated accretion flow (ADAF) around the black hole, and can thus be used to probe the innermost regions 
around the black hole. We calculate the optical depth of the ADAF radiation field to the TeV photons and find it 
strongly depends on the spin of the black hole. We find that transparent radii of 10 TeV photons are of 5/?s and 
13/?s for the maximally rotating and non-rotating black holes, respectively. With the observations, the calculated 
transparent radii strongly suggest the black hole is spinning fast in the galaxy. TeV photons could be used as a 
powerful diagnostic for estimating black hole spins in galaxies in the future. 

Subject headings: black hole physics — galaxies: individual M87 



1. INTRODUCTION 

It is generally accepted that accretion onto supermassive 
black holes (SMBHs) is the major source of their growth during 
the periods of the active nuclei duty cycles (Chokshi & Turner 
1992; Yu & Tremaine 2002; Marconi et al. 2004; Wang et al. 
2008). This process spins up the SMBHs, resulting in that the 
most SMBHs should be rapidly rotating (Bardeen 1972; Thorne 
1974; Volonteri et al. 2005; Wang et al. 2006). The only strong 
evidence for rapidly rotating individual extragalactic SMBHs is 
the gravitationally broaden profile of the iron Ka emission line 
in the X-ray spectrum of MCG-6-30-15 (Wilms et al. 2001; 
Fabian et al. 2002; Brenneman & Reynolds 2006). The general 
lack of observational evidence for gravitationally broaden iron 
Ka lines leaves the issue of black hole spin as one of the most 
elusive questions in astrophysics. 

It is well-known that M87 (at a distance of ^ 16 Mpc) con- 
tains an SMBH with a mass of M. = (3.2 ±0.9) x lO'^M© 
based on observations with Hubble Space Telescope (Harm et 
al. 1994; Macchetto et al. 1997). Thanks to the high spatial 
resolution of these observations, structures in the core of M87 
have been resolved on scales of 20pc at optical wavelengths 
(Harm et al. 1994) and ^ 100/?s in radio wavelengths (Junor 
et al. 1999), where R^, = 2GM,/c^ is the Schwarzschild radius, 
G the gravitational constant and c the light speed. The multi- 
wavelength continuum can be explained by a simple ADAF or 
modified model with outflows under a Schwarzschild hole (Di 
Matteo et al. 2000; 2003). Several parameters in these models 
are strongly degenerated, it remains an open question whether 
the SMBH is rapidly rotating. 

The goal of the present paper is to estimate the SMBH spin 
based on the TeV variability detected in M87 by High Energy 
Stereoscopic System (HESS) (Aharonian et al. 2006). An ori- 
gin from jet has been ruled out (Aharonian et al. 2006). The at- 
tenuation of TeV photons by the ADAF radiation field depends 
on the spectral energy distribution (SED) along the ADAF ra- 
dius. This presents us with the opportunity to tackle this issue. 

2. ADVECTION-DOMINATED ACCRETION FLOW IN M87 



Fig. 1 shows the multiwavelength continuum of M87. From 
radio to X-rays, the SED is dominated by emission from the 
accretion flow around the SMBH (Di Matteo et al. 2003). The 
spectrum has a bolometric luminosity of ^ lO^^'erg s"', giving 
the accretion flow an Eddington ratio of ^ 10"''. This implies 
that an advection-dominated accretion flow is at work in the 
galaxy (Reynolds et al. 1996; Di Matteo et al. 2003). For sim- 
plicity, we use a self-similar solution of the ADAF (Narayan 
& Yi 1994) to determine the temperatures of both protons and 
electrons (Narayan & Yi 1995) in order to calculate the emer- 
gent spectrum. The self-similar solution of the ADAF is char- 
acterized by a factor /, which describes the fraction of the re- 
leased gravitational energy that is advected into the black hole. 
The tendency of / — > 1 has been verified by extensive numerical 
calculations of the global solution of the ADAF (e.g. Narayan 
et al. 1997; Manmoto et al. 1997). We define m = i]Mc^/L^i, 
where M is the accretion rate, ri{= 0.1) the radiative efficiency 
and LEdd the Eddington luminosity. There are five parameters 
{a,j3^m,M,^Rin) in the ADAF model, where a is the viscosity, 
/?in the inner radius of the ADAF, (3 = f'gas/ (^gas + ^fi), ^gas and 
Pb are the gas and magnetic pressures, respectively. We refer 
Schwarzschild and Kerr holes to those with the specific angular 
momentum a = and a = 1, respectively, where a = J /GM\c~^ 
and J is the angular momentum of the hole. Spins are directly 
related to the inner radius R\n , which can be determined by 
fitting observed continuum. For Schwarzschild and Kerr holes, 
the inner radius R\a = 6R0 and Rg, respectively, where Rg = Rs/2 
is the gravitational radius. We begin by fitting the observed 
multiwavelength SED of M87 to determine the best-fit param- 
eters of the accretion flow. Then, we find the optical depth of 
the ADAF radiation to TeV photons. We include synchrotron, 
multiple Compton scattering and free-free emission processes 
in this calculations. We follow the numerical scheme for the 
emergent spectrum described in Manmoto et al. (1997) (see 
also Yuan et al. 2000) in this Letter. 

We used the black hole mass measured by HST in fitting 
the continuum. The viscosity and advection factors of a and 
/ are actually absorbed into one parameter in the self-similar 
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Fig. 1. — The spectral energy distribution of M87 from radio to TeV bands. The three radio data points are from VLBI nuclear flux density measurements 
(Pauliny-Toth et al. 1981; Spencer et al. 1986; Baath et al. 1992 in increasing frequency order). The near IR 10^ nuclear flux measurement is the best IR flux limit 
(a resolution of ~ 0".46) from Gemini observations (Perlman et al. 2001). The optical nuclear continuum measurement is from //ST observations with a resolution 
of 0" . 15 (Harms et al. 1994). The black line with an open circle is the Chandra spectrum (0.2— 10 keV) (Di Matteo et al. 2003). The color lines are different ADAF 
models for Kerr (red) and Schwarzschild (blue) black holes. The black dashed hne is the spectrum from ADAF model with axeV = 0.65 and ni = 1.1 X 10"^. 



solution. Though the viscosity a is highly unknown, only 
their combination of a and / is important for the SED model 
(Narayan & Yi 1995). The final optical depth to TeV photons is 
directly dependent on the radial SED rather than a and /, so we 
take typical values of a = 0.3 and / = 0.99. The ADAF spec- 
trum is characterized by a sharp peak of synchrotron emission 
from moderately relativistic electrons, multiple bumps from in- 
verse Compton scattering and a high energy peak with a cut- 
off originating from bremsstrahlung emission from the very hot 
plasma in the ADAF. The roles of [3 and m are clear in fit- 
tings of the continuum (Narayan & Yi 1995) if the inner ra- 
dius of the ADAF is fixed. For a plasma with equipartition 
between the gas and magnetic fields, /?eq = 0.5. For f3 > /3eq, 
the plasma is weakly magnetized and for /3 < /3eq the plasma 
is strongly magnetized. Usually f3 = 0.5 works generally, but 
it is needed to be elaborately adjusted for specific fittings. For 
a model with a fixed Rin, the synchrotron peak is primarily de- 
termined by M,, while the Compton bumps and X-ray peak are 
sensitive to m and f3. Fig. 1 shows the best fittings ADAF 
models with Rin = l,6Rg. The Kerr model, with (3 = 0.35 and 
m = 5.5 X 10""*, is consistent with the entire spectrum from ra- 
dio to X-rays, while the Schwarzschild model with (3 = 0.35 
and m = 2.0 x 10"^ can match the data at only radio, near IR 
and optical wavelengths, but always not with the X-ray data. 
The differences between the two models originate from the fact 
that the Kerr black hole has higher radiation efficiency than the 
Schwarzschild. For the same total luminosity, the accretion rate 
is lower for the Kerr hole than the Schwarzschild, resulting in a 
thinner medium of the accreting gas and then reduced free-free 
emission in X-ray band. We note, however, that modified mod- 
els of the ADAF might in principle produce a variety of spectra 
(Di Matteo et al. 2003) and that measurements of the SMBH 
spins strongly depend upon m. We present a new method for 



revealing the SMBH spin using the constraints from TeV vari- 
ability. 

Uncertainties in the ADAF fit model results are primarily 
governed by m and for models with fixed a and /. In partic- 
ular, the modeled X-ray spectrum strongly depends on both m 
and Ri„. It is difficult to estimate these uncertainties, however 
the upper limit of Ri„ will be constrained by the upper limit of 
the region responsible for the observed TeV variability. 
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Fig. 2. — The geometric scheme. For simphcity, we assume that the TeV 
photons originate from a point P in the XOY— plane. R is the distance of the 
TeV photons to the hole, Sd is the radius of the ADAF. © is the viewing angle. 

3. OPTICAL DEPTH OF TEV PHOTONS AND HOLE'S SPINS 

We calculate the number density of soft photons from the 
ADAF as a function of position P{R, 9,0). Fig. 2 shows the ge- 
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ometric scheme for the ADAF disk. The vertical structure of the 
ADAF is unknown, but it is generally assumed that the vertical 
density has an exponential profile (e.g. Manmoto et al. 1997). 
We follow this with the assumption that most of the emitted 
photons originate from the mid-plane of the accretion disk. 
With simple manipulation, we have cos(POD) = sin 8 cos 0, and 
the distance 

dl = R^ + Rl-2RRj:)SmQcos4>. (1) 

The cosine of the angle between two interacting photons is then 

R-Rd sin cos 
M= ^ ■ (2) 

We then have the number density of photons from the ADAF, 
given by 



nph(Q,R,e,RB)- 



F(RD,e) 



— I cose, 

do. 



(3) 



where FiRu,e) is the upper stream flux of photons from a radius 
Rd, the factor of (^R/dr,^ cos 8 is the projected area element in 
the direction DP, is the mass of the electron and e is the 
photon energy in units of nieC^. 




Fig. 3. — The SEDs at 4 and 3 different radii for Kerr (red) and 
Schwarzscliild (blue) black holes with the parameters obtained from the con- 
tinuum fitting, respectively. 

Fig. 3 shows the SEDs for both Kerr and Schwarzschild 
black holes at four and three different radii, respectively. The 
SEDs demonstrate that gravitational energy is dissipated in a 
more compact region for a Kerr hole than a Schwarzschild. This 
is the key point that provides an opportunity to probe the black 
hole spin. We stress here that we neglect general relativistic 
(GR) effects in this paper, however the release of gravitational 
energy via viscous dissipation will be concentrated in a more 
compact region if we include the effects. This, in fact, increases 
the robustness of the conclusions presented here. More sophis- 
ticated model including GR effects will be given in a future 
paper 

TeV photons can divulge details of the ambient radiation field 
through their dilution via pair production. Pair production can 
not be avoided if energies of two colliding photons ei and £2 (in 
units of MeC') satisfy the relation 

2 

(l-/l)(l-v2) 

where n is the cosine of the angle between the two photons, and 
V is the velocity of the positrons and electrons in the center of 
momentum frame in the units of c. The optical depth of the ra- 
diation field to TeV photons can be obtained by integrating over 



the radius R along the propagation path of the TeV photons to 
infinity and direction of incident soft photons, i.e.. 



T^JQ,R,e)-- 



dR 



27T 



RodRo / d<j) / de' 



(5) 



a^^(e, e', fi)npt^(Q,R, e',/?D), 

where (T^-y(e, e', /i) is the cross section of the two colliding pho- 
tons (Gould & Schreder 1967), e, is the threshold energy of the 
interacting photons and Rom is the outer radius of the ADAF. It 
should be noted that t~^^ strongly depends on the incident angle 
between the two photons. 

The dependence of t^~^ on Rin (and thus the SMBH spin) 
can be qualitatively predicted. For a given bolometric luminos- 
ity, most of the gravitational energy dissipated during accretion 
onto a rapidly black hole is released in a more compact region 
than in a system containing a Schwarzschild black hole. This 
leads to a much steeper r^^ - R relation around the Kerr black 
hole horizon and facilitates the escape of TeV photons from 
considerable smaller radii. It is thus expected that the timing 
information from TeV photons can be used as a diagnostic of 
the spatial distributions of the soft photon field surrounding the 
SMBH and thus its spin. 




Fig. 4. — The optical depth of the ADAF radiation field to 10 TeV photons. 
The solid and dashed lines are for © = 0° and 30°, respectively. The radio jet 
has an orientation angle of 30° with respect to an observer's sight line, so we 
assume the ADAF is viewed at = 30°. We only plot t^^. to lOTeV photons, 
because Tj-f < 1 for most radii for photons < ITeV. 



Fig. 4 shows the optical depth of the radiation field to 
10 TeV photons for the best fitting ADAF models for both 
Schwarzschild and Kerr black holes. We fix Rqm = lO^Rg- We 
note that the final results are not sensitive to 7?ouf t^-j is sensi- 
tive to radii, but not to the viewing angle. As we argue quali- 
tatively. Fig. 4 shows that the t^^-R relation is much steeper 
for Kerr holes than the Schwarzschild. For TeV photons to be 
able to escape from the radiation field of the ADAF, we require 
< 1. We define the transparent radius, R^, as the radius 



at which t^^ = 1 



1 IRg and 26Rg for Kerr and 



We find R. 

Schwarzschild black holes, respectively. The TeV variability 
timescale in M87, observed by HESS, is of ~ 2days, imply 
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ing a compact emission region of /Jxev < cVAtjcv ~ 10 VR^, 
where V =l/T (l-VcosO/c) is the Doppler factor and F = 

2-1-1/2 

l-(y/c) the Lorentz factor of a jet with a velocity of 

If the jet is viewed at 9 = 30° (Bicknell & Begelman 1996; 
Aharonian et al. 2006), the maximum Doppler factor is ©max = 
2. Even if the TeV photons are emitted from the base of the rel- 
ativistic jet, we have /?TeV < 20 2?2/?g, where I>2 = I'max/2. On 
the necessary condition of /?TeV > Rc, the HESS observations 
suggest the presence of a rapidly rotating black hole in M87 
in term of Rc = 26/?g > /?TeV for a Schwarzschild hole. Other- 
wise, the observed TeV region would be optical thick to lOTeV. 
Increasing the black hole spin, a, the transparent radius R^ de- 
creases. In order for the transparent radius to equal the TeV 
emission radius, we require a black hole spin of axev = 0.65 for 
m = 1 . 1 X 1 0~^ . We note that the X-ray spectrum can not be well 
fit with a = . 65 and w = 1 . 1 X 1 0~^ , however we show the SED 
in Figs 1 and 4 to demonstrate how sensitive the transparent ra- 
dius and the SED are to black hole spin. Fig 4 shows a region 
between 1 > a > axeV. in which the exact spins should be, most 
likely have the maximum spin if the TeV region has a radius of 
lOR,. 

We note that although the t^^ - R relation becomes steeper 
with spin, R^ is in fact insensitive to spins of a < axeV, making 
this technique less useful for estimation of their spins in slowly 
rotating black holes. The absence of any evidence for a poten- 
tial break between 10-20TeV means Tj~^ <C 1 for these TeV 
photons, namely the spin a > axeV- Otherwise, atR= lORg for 
a Schwarzschild hole, r-y^ = 5 ~ 6 and the TeV photons will be 
attenuated by a factor of e'^"''' w 150 ~ 400. This strongly con- 
flicts with results from the TeV observations. The conclusion 
does, however, depend on the fitting of the multiwavelength 
continuum, which could be contaminated by other sources, for 
example, the inner jet. The current data and model do not allow 
us to constrain the spin to better than a > 0.65, however, fu- 
ture improved multiwavelength data and models will produce 
tighter constraints on how fast the spin is. 

A rapid rotation is generally consistent with the statistic re- 
sults from the Sloan Digital Sky Survey (Wang et al. 2006) 
and with justification from radio emission in term of Blandford- 
Znajek effect (Newman et al. 2007). Since SMBH growth is 
primarily accretion, the spin angular momentum of the black 
hole is Ukely to be dominated by its accretion history. If the 
accretion is random, however, as suggested by King & Pringle 
(2007), the net spins of SMBHs should tend to zero. The spin 
evolution of SMBHs will thus uncover the accretion modes (co- 
rotating, retrograde or random accretion) in their history from 



measuring spins of local SMBHs. 

4. CONCLUSIONS AND DISCUSSIONS 

The TeV variabihty as a new independent ingredient allows 
us to disentangle the inner radius from other parameters in mod- 
eling the multiwavelength continuum. We fit the multiwave- 
length continuum of M87 and use the best-fit parameters of 
the ADAF to calculate the optical depth of the ADAF radiation 
field to TeV photons. This optical depth is sensitive to the spin 
of SMBHs and combined with the TeV variability constraints, 
we conclude that the SMBH must be rapidly rotating in M87. 
Whatever mechanism produces the TeV photons in M87, these 
energetic photons provide a potentially powerful probe of the 
SMBH spin. 

We should note some uncertainties of the present results, 
which fully depends on the theoretical model of ADAF, and 
contamination of the radiation in the central region from its host 
and jet. The influence of potential outflows from the ADAF 
should be examined in a future paper. Although we neglect the 
GR effects and the global structure of the ADAF in this paper 
for simplicity, the addition of the effects will make the conclu- 
sion of this work robust, because the gravitational energy is dis- 
sipated in more compact regions due to the deeper gravitational 
potential than in the Newtonian case. The resulting t^-^ -R re- 
lation will then be steeper, reducing the transparent radius from 
which TeV photons can escape. Accurate measurements of the 
black hole spin rely on global calculations of the ADAF with 
general relativity and the timescale measurements of TeV vari- 
ability to confine the dimensions of the TeV region. Future ob- 
servations of other sources with HESS may reveal similar TeV 
variability to M87, particularly in low luminosity active galac- 
tic nuclei (LLAGNs), such as radio galaxies and low ionization 
nuclear emission regions (LINERs). This would lead to further 
estimates of the spins of SMBHs in nearby galaxies. Addition- 
ally, the present results suggest that any mechanism responsible 
for production of the TeV photons should be located outside a 
region with r^-y > 1, constraining on the radiation mechanism 
of TeV photons. 
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